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Abstract. The chemistry of CH 3 0H and H 2 CO in thir- 
teen regions of massive star formation is studied through 
single-dish and interferometer line observations at sub- 
millimeter wavelengths. Single-dish spectra at 241 and 
338 GHz indicate that T rot = 30 - 200 K for CH 3 OH, 
but only 60 - 90 K for H 2 CO. The tight correlation be- 
tween T rot (CH 3 OH) and T ox (C 2 H 2 ) from infrared absorp- 
tion suggests a common origin of these species, presum- 
ably outgassing of icy grain mantles. The CH3OH line 
widths are 3 — 5 km s _1 , consistent with those found ear- 
lier for C 17 and C 34 S, except in GL 7009S and IRAS 
20126, whose line shapes reveal CH3OH in the outflows. 
This difference suggests that for low-luminosity objects, 
desorption of CHaOH-rich ice mantles is dominated by 
shocks, while radiation is more important around massive 
stars. 

The wealth of CH3OH and H 2 CO lines covering a large 
range of excitation conditions allows us to calculate radial 
abundance profiles, using the physical structures of the 
sources derived earlier from submillimeter continuum and 
CS line data. The data indicate three types of abundance 
profiles: flat profiles at CH 3 OH/H 2 ~ 10" 9 for the cold- 
est sources, profiles with a jump in its abundance from 
~ 10~ 9 to ~ 10~ 7 for the warmer sources, and flat pro- 
files at CH 3 OH/H 2 - few 10~ 8 for the hot cores. The 
models are consistent with the f» 3" size of the CH3OH 
107 GHz emission measured interferometrically. The lo- 
cation of the jump at T « 100 K suggests that it is due 
to evaporation of grain mantles, followed by destruction 
in gas-phase reactions in the hot core stage. In contrast, 
the H 2 CO data can be well fit with a constant abundance 
of a few xlO" 9 throughout the envelope, providing lim- 
its on its grain surface formation. These results indicate 
that T rot (CH3OH) can be used as evolutionary indica- 
tor during the embedded phase of massive star formation, 
independent of source optical depth or orientation. 

Model calculations of gas-grain chemistry show that 
CO is primarily reduced (into CH 3 OH) at densities 
riH~10 4 cm -3 , and primarily oxidized (into C0 2 ) at 



higher densities. A temperature of w 15 K is required 
to keep sufficient CO and H on the grain surface, but 
reactions may continue at higher temperatures if H and 
O atoms can be trapped inside the ice layer. Assum- 
ing grain surface chemistry running at the accretion rate 
of CO, the observed abundances of solid CO, C0 2 and 
CH3OH constrain the density in the pre-protostellar phase 
to be rtn~ a few 10 4 cm" 3 , and the time spent in this 
phase to be < 10 5 yr. Ultraviolet photolysis and radiolysis 
by cosmic rays appear less efficient ice processing mech- 
anisms in embedded regions; radiolysis also overproduces 
HCOOH and CH 4 . 

Key words: ISM: molecules - Molecular processes - 
Stars: Circumstellar matter; Stars: formation 



1. Introduction 

One of the early stages of massive star forma- 
tion is the "hot core" pha se (Kurtz et al. 2000, 



Macdonald fc Thompson 2000 ) , which is characterized by 
masses of ~ 100 M Q of molecular gas at temperatures 
of > 100 K, leading to a rich line spectrum at submil- 
limeter wavelengths like that of the prototypical Orion 
hot core. The chemical composition of hot cores is quite 
distinct for its high abundances of fully hydrogenated, 
large carbon-bearing molecules (Ohishi 1997), as opposed 
to dark clouds which show predominantly unsaturated 
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species and molecular ions (van Dishoeck & Blake 1998). 
Saturated carbon chains are not expected in large 
quantities by steady-state gas-phase models, which led 
|Charnlcy et al. (1992)| and paselli et al. (1993)| to pro- 
pose models where these molecules are made on dust 
grains, and evaporate into the gas phase when the newly- 
formed star heats up its surroundings. This leads to a 
short period (~ 10 4 yr) when the envelope of the young 
star is rich in complex organic molecules. Grain sur- 
face reactions at low temperature lead to hydrogenation 
if atomic H is able to qu antum tunnel through a ctiva- 
tion barriers as proposed by Ticlcns fc Hagen (1982) . This 
model makes specific predictions for the abundance ra- 
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Spectra of the CH3OH J = 7k — > §k transition at 338 GHz, obtained with the JCMT. 



tios of CO, H2CO, CH3OH and their deuterated versions 
which for a reasonable choice of parameters fit the abun- 
dances observed in the Compact Ridge in Orion as shown 



by Charnley et al. (1997). However, it is unknown if other 
sources can be fit as well. Alternatively, processing by ul- 
traviolet light or energetic particles may change the ice 
composition in the vicinity of young stars. 

One key molecule to test models of hot core chemistry 
is methanol, CH3OH. At temperatures < 100 K, produc- 
tion of CH3OH in the gas phase goes mainly via radiative 
association of CH^ and H2O, an inefficient process yield- 
ing abundances of only ~ 10 -11 relative to H2 (Lee et al. 
1996). Observations of CH3OH toward dark and translu- 
cent clouds by Turner (1998) and Takakuwa et al. (1998) 
yield abundances of ~ 10~ 9 , which suggests that grain 
chemistry operates. Toward the Orion hot core and com- 
pact ridge, CH3OH abundances of ~ 10 -7 have been 



observed by |Blake et al. (1987)| and [Sutton et al. (1995) 
and even higher abundances, ~ 10~ °, have been 
inferred from obs ervations of CH3OH masers (e.g., 
Menten et al. 1986 ). Methanol has also been observed in 



the solid state at infrared wavelengths from the ground 
by, e.g., |Grim et al. (1991)| and [Allamandola et al. (1992) 
and recently by Dartois et al. (1999b), at abundances up 
to 10~ 6 . 

To study the organic chemistry in warm (30 — 200 K) 
and dense (10 4 — 10 7 cm~ 3 ) circumstellar envelopes, this 
paper discusses observations of CH3OH and H2CO lines 
at submillimeter wavelengths. Although gas-grain inter- 
actions are clearly important for methanol, it is unknown 
how efficiently (if at all) surface reactions modify the com- 
position of the ices. In addition, it is still open what 
process returns the molecules to the gas phase: thermal 
heating, shocks or both? Finally, the ice layer must have 
formed at temperatures and densities much lower than 
the present situation, and hence provides a unique fossil 
record of the conditions in the molecular cloud prior to 
star formation. 

The sources are thirteen massive (L = 10 3 — 10 5 L©) 
stars, which are at an early stage of evolution and still 
embedded in 10 2 — 10 3 M© of dust and molecular gas. 
A few of the sources are hot cores by the definition of 
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Fig. 2. Line profiles of CH 3 OH 5_ -» 4_ A+ toward 
GL 7009S and IRAS 20f26. The outflow wing is most 
prominent in GL 7009S. The feature at —25 km s _1 in 
IRAS 20126 is the K = 4 A line. 



Fig. 3. Rotation temperatures of CH3OH (left) and H2CO 
(right) versus excitation temperatures of C2H2 measured 
in infrared absorption. Crosses in the bottom right denote 
typical error margins. 



Kurtz et ah (2000), but most are in an even younger phase 
where most of the envelope is still at low temperatures. 
The physical structure of the soures has been studied by 
van der Tak et ah (2000), who developed detailed temper- 
ature and density profiles. The wealth of CH3OH and 
H2CO lines, combined with the physical models, allows 
the determination of abundance profiles for H2CO and 
CH3OH, which demonstrate that CH3OH evaporates off 
dust particles in the envelope ona~ 10 4-5 yr time scale, 
while H2CO is predominantly formed in the gas phase. 
The excitation and abundance of CH3OH are therefore 
useful evolutionary indicators during the embedded stage 
of star formation. Moreover, by comparing the observed 
amounts of evaporated CH3OH and CO2 to a model of 
grain surface chemistry, we set limits on the density and 
the duration of the pre-protostellar phase. 

2. Observations 

2.1. Single-dish observations 

Spectroscopy of the J = 5 — » 4 and 7^6 bands of 
CH3OH near 242 and 338 GHz was performed with the 
15-m James Clerk Maxwell Telescope (JCMTf] on Mauna 
Kea, Hawaii in May and October of 1995. The beam size 
(FWHM) and main beam efficiency of the antenna were 
18" and 69% at 242 GHz and 14" and 58% at 338 GHz. 
The frontcnds were the receivers A2 and B3i; the backend 
was the Digital Autocorrelation Spectrometer, covering 

1 The James Clerk Maxwell Telescope is operated by the 
Joint Astronomy Centre, on behalf of the Particle Physics 
and Astronomy Research Council of the United Kingdom, the 
Netherlands Organization for Scientific Research and the Na- 
tional Research Council of Canada. 



500 MHz instantaneous bandwidth. Pointing was checked 
every 2 hours during the observing and was always found 
to be within 5". To subtract the atmospheric and instru- 
mental background, a reference position 180" East was 
observed. Integration times are 30-40 minutes for each 
frequency setting, resulting in rms noise levels in T m b per 
625 kHz channel of » 30 mK at 242 GHz to w 50 mK at 
338 GHz. Although the absolute calibration is only cor- 
rect to « 30%, the relative strength of lines within either 
frequency setting is much more accurate. 

In addition to these CH3OH lines, this pa- 
per discusses observations of lines of H2CO, ob- 
tained in a similar manner, which were presented in 
yan der Tak et al. (2000) . We will also use data on both 



molecules for W 3 IRS5 and W 3 (H2O) from the survey by 
Hclmich fc van Dishocck (1997)|, and CH 3 O H J = 7 -> 6 



data from Thompson & MacDonald (1999) 



2.2. Interferometric observations 

Maps of the J K = 3i -> 4 A+ (107013.852 MHz; 
E u = 28.01 K) and ll_i -> 10_ 2 E (104300.461 MHz; 
E u = 157.24 K) lines of CH 3 OH were obtained for four 
sources in 1995-1999 with the millimeter interferometer 
of the Owens Valley Radio Observatory (OVRO)g. This 
instrument consists of six 10.4 m antennas on North- 
South and East- West baselines, and a detailed technical 



description is given in Padin et al. (1991). For the sources 
GL 2591, NGC 7538 IRS1 and NGC 7538 IRS9, data were 
taken in the compact and extended configurations, while 



2 The Owens Valley Millimeter Array is operated by the Cal- 
ifornia Institute of Technology under funding from the U.S. 
National Science Foundation (AST-9981546). 
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Table 1. Fluxes / T MB dV (K km s" 1 ) and FWHM widths (km s" 1 ) of lines observed with the JCMT. 



Line Frequency GL 490 W 33A GL GL S 140 NGC 7538 NGC 6334 GL IRAS W 28 
(K,A/E) a (MHz) 2136 2591 IRS1 IRS1 IRS9 IRS1 7009S 20126 A2 d 
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Line Width 



3.4 5.2 3.0 3.5 2.7 3.9 3.6 5.7 6.9 7.5 6.2 

±0.7 ±1.1 ±1.3 ±0.4 ±0.2 ±0.6 ±0.7 ±0.5 ±1.3 ±1.9 ±0.9 



a When no superscript is given for A-type methanol, the A + and A lines are blended. 
b Blend with K=6 E at 338404.5 MHz, assumed to equal K=-6 E. 
c Possible contribution from K=l A~ at 341415.6 MHz. 

d This source is also known as G 5.89-0.39; 338.5 GHz data are from Thompson & Macdonald 1999. 



for the Southern source W 33A, a hybrid configuration 
with long North-South but short East- West baselines was 
also used to improve the beam shape. Antenna spacings 
range from the shadowing limit out to 90 kA, correspond- 
ing to spatial frequencies of 0.018 — 0.44 arcsec -1 . The 
observations were carried out simultaneously with contin- 
uum observations, which have been presented in van der 
Tak et al. (1999, 2000). We refer the reader to these papers 
for further observational details. 



3. Results 

3.1. Submillimeter spectroscopic results 

Figure |l| presents the calibrated JCMT spectra, reduced 
using the IRAM CLASS package. All lines are due to ei- 
ther A- or E-type CH3OH, while some features are blends 
of lines. After subtraction of a linear baseline, Gaussian 
profiles were fitted to extract line parameters. The free 
parameters in these fits were the line flux J T m \^dV , the 
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Table 2. Beam-averaged column densities and excitation temperatures of CH3OH and H2CO, and abundances inferred 
for a power law envelope using spherical Monte Carlo models. 



Source 


iV(CH 3 OH) 
cm -2 


JV(H 2 CO) 
cm -2 


T(CH 3 OH) 


T(H 2 CO) 


CH3OH/H2 


H 2 CO/H 2 (c) 




K 


K 


10" 9 


10" 9 


W 3 IRS5 


1.6(14) 


7.8(13) 


73 


78 


0.4 


3 


GL 490 


3.6(14) 


4.3(13) 


24 


94 


1.0 


1 


W 33A 


2.0(15) 


1.2(14) 


155 


88 


3.1 / 90 (6) 


4 


GL 2136 


4.4(14) 


4.5(13) 


143 


76 


0.9 


8 


GL 7009S (a) 


3.9(15) 




8 




0.7 




GL 2591 


1.2(15) 


8.0(13) 


163 


89 


2.6 / 80 (6) 


4 


S 140 IRS1 


2.5(14) 


8.0(13) 


41 


60 


1.2 


5 


NGC 7538 IRS1 


2.2(15) 


1.9(14) 


189 


87 


2.0 / 60 (6) 


10 


NGC 7538 IRS9 


6.5(14) 


1.0(14) 


29 


82 


2.3 


10 


W 3 (H 2 0) 


7.5(15) 


5.0(14) 


203 


181 


5.9 


3 


NGC 6334 IRS1 


3.8(16) 


1.6(15) 


213 


193 


24. 


7 


IRAS 20126 (a) 


2.2(15) 




139 




2.6 




W 28 A2 


4.2(15) 




43 




12. 





(a) : Only CH3OH J = 5 4 observed. 

(b) : Abundances in cold (T < 90 K) and in warm (T > 90 K) gas. 

(c) : From van der Tak et al. 2000. 



FWHM line width and the central velocity. The retrieved 
line fluxes and widths are listed in Table [l], while the cen- 
tral velocities are consistent with those found for C 17 
and C 34 S in |van der Tak et al. (2000)] 



For most sources, the line widths are also consistent 
with those found for C 17 and C 34 S. However, for the 
two lowest-luminosity sources, GL 7009S and IRAS 20126, 
the methanol lines are much broader than those of C 17 
and C 34 S. The CH3OH line profiles toward these sources 
(Fig. |^) suggest an origin in their molecular outflows, espe- 
cially in GL 7009S. The metha nol emission toward y oung 
low-mass stars such a s L 1157 (Bachillcr et al. 1995 ) and 
NGC 1333 IRAS 4A flBlakc et al. 1995| ) is also dominated 
by their outflows. The other sources studied here pos- 
sess outflows, but those do not seem to be important 
for the CH3OH emission. As we will show below, the ob- 
served gas-phase CH3OH is mostly due to evaporation of 
grain mantles. The CH3OH line profiles thus suggest that 
around low-mass stars, the desorption of icy grain mantles 
is dominated by shocks, while radiative heating dominates 
for high-mass stars. 

The sources exhibit clear variations in their methanol 
spectra, both in absolute line strength and in the relative 
strength between lines. Variations in absolute line strength 
can be due to differences in molecular abundances, the size 
of the emitting region or distance, while the line ratios 
reflect differences in the temperature and density of the 
emitting gas. Before modeling the observations in detail, 
we perform a quick-look examination using rotation dia- 
grams, as in Blake et al. (1987) and Hclmich et al. (1994) 



be described by a single temperature, the "rotational tem- 
perature" . Since the envelopes of these stars are known 
to have temperature and density gradients, the rotation 
temperature gives information where the molecule is pref- 
erentially located in the envelope. The nondetection of 
13 CH30H lines limits the optical depth of the lines to < 3, 
except in the case of NGC 6334 IRS1, where the strongest 
lines are probably saturated. 

Table || presents the results of the rotation diagram 
analysis, which are averages over the JCMT beam. For- 
mal errors from the least-square fit to the data and cali- 
bration errors combine into uncertainties of 20% in both 
parameters. The results for IRAS 20126, W 3 (H 2 0) and 
NGC 6334 IRS1 are more uncertain, to ps 50%, because 
the derived T rot is close to the highest observed energy 
level. The OVRO maps presented in § 3.2 indicate that 



This analysis assumes that the emission is optically thin 
and that the population of the molecular energy levels can 



the CH3OH column densities are beam diluted by a factor 
of ~ 25; this information is not available for H2CO. Al- 
though the column densities of the two species are roughly 
correlated, with AT(CH 3 OH) - 10 x A^(H 2 CO), their exci- 
tation temperatures behave differently. While T ro t(H 2 CO) 
is rather constant from source to source, 60 — 90 K, except 
toward the hot cores W 3 (H 2 0) and NGC 6334 IRS1, 
the CH3OH excitation temperatures span the full range 
10 - 200 K. 

The spectral line survey of 

Thompson fc MacDonald (1999)| toward W 28 A2 
includes the CH3OH J = 7 — > 6 band, and their analysis 
gives T rot = 49 K and N = 3.8 x 10 15 cur 2 , in good 
agreement with our results from the J = 5 — » 4 band. Our 
values of N(R 2 GO) toward W 3 (H 2 0) and NGC 6334 
IRS1 are a factor of ~10 higher than those found by 
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Fig. 4. Maps and spectra of CH3OH line emission, made with the OVRO array. Contour levels for GL 2591 (top row) 
are: 60, 180, 300 mJy beam" 1 for both images; for NGC 7538 IRS1 (middle row): 300, 900, 1500 mJy beam" 1 (left) 
and 400, 2400 and 4400 mJy beam" 1 (right); and for W 33A (bottom row): 100, 200 mJy beam" 1 (left) and 100, 300, 
500, 700 mJy beam -1 (right). In the spectrum of NGC 7538 IRS1, the brightness of the — 10o E line has been 
increased by a factor of 4. 



Mangum fc Woottcn (1993)| , probably due to the smaller 
JCMT beam size. 

Figure [?] compares the rotational temperatures of 
CH3OH and H2CO to the C2H2 excitation temperatures 
measured by Lahuis fc van Dishoeck (2000)| in absorption 
at 13.7 /im for the same lines of sight. Lacking a dipole mo- 
ment, C2H2 is a clean probe of temperature. The tight cor- 
relation between the CH 3 OH and C2H2 excitation temper- 
atures implies that the two species trace the same gas, and 
that source orientation and optical depth effects do not in- 
fluence the infrared data. The only exception is W 33A: 
this massive source becomes optically thick at 14 /im be- 
fore the warm gas is reached. In contrast, the H2CO tem- 
perature is not correlated with that of C2H2: omitting 
W 33A, the correlation coefficient is 7.8%, versus 98.3% 
for CH3OH. This difference indicates that in these sources, 
H2CO is not chemically related to CH3OH and C2H2, un- 
like in low-mass objec ts, as the high H 2 CO abundance 
in the L 1157 outflow flBachillcr fc Perez Gutierrez 1997|) 



and the detection of HDCO and D 2 CO in IRAS 16293 



(van Dishoeck et al. 1995; Ceccarelli et al. 1998) indicate 



3.2. Interferometric maps of CH3OH line emission 

Figure ^ presents the OVRO observations, as maps of 
the emission integrated over the line profile. These maps 
were obtained by a gridding and fast Fourier transform of 
the (u, w)-data with natural weighting, and deconvolution 
with the CLEAN algorithm. The self-calibration solutions 
obtained for the continuum data were used to suppress 
phase noise introduced by the atmosphere. Also shown 
are spectra at the image maxima, extracted from image 
cubes at full spectral resolution and coverage. In these 
maps and spectra, the continuum emission has not been 
subtracted. Only in the case of NGC 7538 IRS1, the con- 
tinuum brightness is comparable to the line brightness; for 
the other sources, the continuum can be neglected. 
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Source 


Line 






LSR Velocity 
km s _1 


AV (FWHM) 
km s _1 


Peak T B 
K 


Source Size 
arcsec 


Beam Size 
arcsec 


GL 2591 


3i -> 


4 A+ 




-5.8 


3.5 


4.3 


3.4 x 2.7 


3.4 x 2.9 




ll-i 


-> 10-2 


E 


-5.7 


4.0 


2.8 


3.1 x 2.3 


3.8 x 3.0 


W 33A 


3i -> 


4 A+ 




39.5 


3.5 


23.3 


3.9 x 3.0 


5.6 x 3.4 




ll-i 


-> 10-2 


E 


39.1 


4.3 


4.4 


6.4 x 2.2 


5.7 x 3.5 


NGC 7538 IRS1 


3i 


4 A+ 




-60.5 


1.4 


63.6 




3.9 x 3.2 










-58.6 


1.8 


108. 




3.9 x 3.2 










-56.1 


1.9 


164. 




3.9 x 3.2 










-53.9 


2.5 


29.2 




3.9 x 3.2 




ll-i 


-> 10-2 


E 


-59.0 


2.8 


4.2 


1.8 x 0.7 


4.0 x 3.3 



The maps show single, compact emitting sources, 
which coincide within the errors with their counterparts 
in 86 — 230 GHz continuum emission (van der Tak et al. 
1999, 2000). No CH3OH line emission was detected toward 
NGC 7538 IRS9 to the la noise level of w 1 K. In the 
case of W 33A, a binary source at these wavelengths, the 
CH3OH emission is associated with the source MM1 from 



van der Tak et al. (2000) and with the infrared source. Al- 
though the CH3OH emission is seen to be slightly ex- 
tended toward the West, no emission was detected at the 
position of the other continuum source, MM2. 

Table ^| lists various parameters of the lines observed 
with OVRO. Deconvolved source sizes were obtained by 
fitting two-dimensional Gaussian profiles to the integrated 
emission maps; other parameters by fitting Gaussians to 
the presented spectra. The continuum emission was sub- 
tracted from the data before this fitting procedure. 

The measured central velocities and line widths are 
consistent with the values measured with the JCMT. To 
check if the two telescopes indeed trace the same gas, 
we have calculated the line strengths expected from the 
JCMT rotation diagrams. For the sources where CH3OH 
is detected with OVRO, we predict a brightness of 15 — 
30 K in the CH 3 OH 107 GHz line and of 10 - 15 K in 
the 104 GHz line, in good agreement with the data. For 
NGC 7538 IRS9, which has a low CH 3 OH rotation temper- 
ature, we expect Tb ~ 60 K in the 107 GHz line. The lack 
of CH3OH emission in the OVRO data suggests that there 
is a gradient in T ex from k 30 K on a 15" scale (traced by 
the JCMT but resolved out by OVRO) to ps 100 K within 
the 3" OVRO synthesized beam. For the other sources, the 
bulk of the envelope may already be heated to w 100 K 
and up. 

3.3. Maser emission 

While the Vlsr = —53.9 km s^ 1 component in the 
107 GHz line profile observed toward NGC 7538 IRS1 
with OVRO has a brightness and a width consistent 



with the JCMT data of this source, the brightness of 
the other velocity components is a factor of 2 — 5 higher 
than expected from the JCMT data. This difference sug- 
gests that the observed emission is amplified by a maser 
process. This possibility is supported by the line pro- 
file with three narrow velocity peaks (Figure [|). Maser 
emission has been observed toward this source in several 
lower-frequency methanol lines, by Wilson et al. (1984) 



and Batrla et al. (1987), at velocities consistent with the 
two strongest components reported here. Using VLBI, 
Minicr et al. (1998) resolved the maser emission in the 
5i ^ 6 A+ (6.67 GHz) and 2 -> 3_i E (12.2 GHz) 
lines into eleven spots, with sizes of < 10~ 3 arc seconds 



each. The spots at velocities —57 



-62 km s lie in 



an outflow, but the kinematics of brightest masers, those 
at ps —56 km s _1 , are consistent with an origin in a disk. 
These two lines as well as the 107 GHz line are transitions 
to the K— ground state (so-called Class II masers), which 
r equire a strong co ntinuum radiation field to be pumped 
( ICragg et al. 1992|) . In the case of NGC 7538 IRS1, this ra- 
diation is provided by the H II region, consistent with the 
fact that among the sources studied here with OVRO, only 
NGC 7538 IRS1 displays maser emission, which is also the 
source with the stro ngest H II region. For the lm — > In 



line of H2CO, where Rots et al. (1981) observed maser ac 



tion at similar velocities, t he same pump mechanism op- 
erates, as demonstrated by 



Boland fc de Jong (1981/ 



The brightness ratio of the 6.67 
methanol maser lines measured by 
is w 1 0. Comparing this number 
els by ISobolcv et al. (1997) 
10 16 



to 12.2 GHz 



Minier et al. (1998) 



wc 



with the mod- 
find Af(CH 3 OH)~ 



cm 



Boland fc de Jong (1981 



This value is similar to A^^CO) from 
which is a lower limit because 



Rots et al. (1981) did not resolve the emission spatially. 



The ratio Af(H 2 CO)/Af(CH 3 OH) for the disk of NGC 7538 
IRS1 is therefore > 1, which is significantly larger than in 
the extended envelopes. 
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4. Abundance profiles of CH 3 OH and H 2 CO 

In this section, we proceed to interpret the CH 3 OH and 
H2CO observatio ns in terms of the detailed physical mod- 
els developed by van der Tak et al. (2000) . These models 
use a power law density structure, n = no(r/ro)~ a . Spe- 
cific values for the density gradient a, the size scale rn 
and the density scale rin for every source can be found 
in van der Tak et al. (2000) . These values, as well as tem- 



perature profiles, have been derived from dust continuum 
and CS, C 34 S and C 17 line observations at submillimeter 
wavelengths. 

4-1- Constant- abundance models 

Using the density and temperature profiles, we have 
first attempted to model the data using constant abun- 
dances of H2CO and CH3OH throughout the envelopes. 
This assumption is motivated by the fact that in mod- 
els of pure gas-phase chemistry, the abundances of 
H2CO and CH3OH do not change much in the range 
20 — 100 K. Molecular excitation and radiative trans- 
fer are solved simultaneously with a computer pro- 
gram based on the Monte Carlo method, written by 
Hogerheijde fc van der Tak (2000) The emission in all 
observed lines is calculated, integrated over velocity, and 
convolved to the appropriate telescope beam. Comparison 
to observations proceeds with the x 2 statistic described in 



van der Tak et al. (2000) 



The collisional rate coefficients for CH3OH have been 
provided by M. Walmsley (1999, priv. comm.), and 



are based on the experiments by [Lees fc Hague (1974) 



A detailed de scription of the coefficients is given by 
Turner (1998) . Following Johnston ct al. (1992) , we have 
set the AK = 3 rates to 10% of the AK = 1 values. 

The results of the model calculations are shown in 
Fig. H an d the derived best-fit abundances are given 
in column 6 of Table |[ In the figure, the observed 
and modeled line fluxes have been converted to column 



densities following Hclmich ct al. (1994) . Three types of 
sources can be distinguished. Most sources can be fitted 
with CH30H/H2^ 10~ 9 , similar to the values found for 
dark and translucent clouds (Turner 1998, Takakuwa et 
al. 1998), while the "hot core"-type sources W 3(H20), 
NGC 6334 IRS1 and W 28 A2 require abundances of a few 
times 10~ 8 . However, for three sources, no single methanol 
abundance gives a good fit. 

4-2. Jump models 

For the sources W 33A, GL 2591 and NGC 7538 IRS1, 
the lines from energy levels > 100 K above ground require 
significantly higher abundances than the lower-excitation 
lines. These results suggest that the abundance of CH 3 OH 
is higher in the warm, dense gas close to the star than in 
the more extended, cold and tenuous gas. From the loca- 
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Fig. 5. Results of Monte Carlo models for CH3OH, with 
'cold' sources at the top, 'jump' sources in the middle and 
'hot cores' at the bottom. Filled symbols are line strengths 
predicted by the best-fit constant-abundance models (cir- 
cles) and jump models (squares) for W 33A, GL 2591 and 
NGC 7538 IRS1. Open symbols are JCMT data, with cir- 
cles indicating detections and triangles upper limits. The 
data on W 28 A2 include those of Thompson & Macdonald 
1999. 
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tion of the break in the abundance profile, it seems likely 
that evaporation of grain mantles, which also occurs at 
~ 100 K, plays a role. As a simple test of ice evapora- 
tion, we have run models for these three sources where 
the CH3OH abundance follows a step function. In these 
"jump" models, the CH3OH abundance is at a low level, 
the "base level", far from the star, while at a threshold 
temperature, the abundance surges to a high value, the 
"top level". The situation is sketched in Figure ^. The 
motivation for this model is that if methanol is present 
in icy grain mantles, its abundance will increase strongly 
when these ices evaporate. For the temperature threshold, 
we take 90 K, which is where water ice, the most refrac- 
tory and most abundant component of the grain mantles, 
evaporates in ~ 10 yr ( Sandford fc Allamandola 1993 ). 

Alternatively, the ice may be desorbed in a weak shock 
associated with the molecular outfl ows these sources are 
known to have: Jones et al. (1996) calculate that a local 



grain-grain velocity dispersion of « 2 km s _1 is sufficient 
to shatter ice material. Shock desorption of methanol is 
known to be important in low-mass protostars, for exam- 



ple L 1157 flBachiller et al. 1995j ) and NGC 1333 IRAS 4A 
(Blake et al. 1995). The large line widths and line profiles 
measured here for the low-luminosity sources GL 7009S 
and IRAS 20126 indicate that shock desorption plays a 
role there; for the other sources, thermal effects probably 
dominate. Since CH3OH is confined to the refractory (po- 
lar) component of the ice, it is not necessary to consider 
the effect of a slowly rising temperature in the hot core 



region, as Viti fc Williams (1999) did 



For the base level, we take the abundances found in 
the previous section, which were constrained mostly by 
the low-JsT lines with E < 100 K. We have considered 
jump factors of 3, 10, 30, 100 and 300 without further 
iteration. Between these models, jumps by factors of ~ 
30 give the best match to the high-excitation lines. The 
results of these models are plotted as squares in Fig. [| and 
listed in column 6 of Table @. 

Further constraints on the abundance profile of 
CH3OH in our sources may be obtained by comparing 
the models to the OVRO data. The points in Fig. |^ are 
the OVRO visibility data of W 33A and GL 2591. The ob- 
servations of the 107 GHz line have been integrated over 
velocity and binned in annuli around the source position. 
Superposed are model points for the constant-abundance 
model and for the jump model. The data do not particu- 
larly favour one model or the other, perhaps because this 
transition traces mostly cold gas. Collisional rate coef- 
ficients up to J = 11 are eagerly awaited, so that the 
104 GHz line can be modeled as well. 

5. Discussion 

5.1. Relation to solid-state observations 

The enhancement of the methanol abundance at high 
temperatures is very likely related to the evaporation 




100 150 
Temperature (K) 



Fig. 6. Density and temperature structure (solid line) and 
abundance profiles (dotted and dashed lines) of the two 
models for CH3OH in the source GL 2591. 



of solid methanol, which has been observed at in- 
frared wavelengths toward many of the sources stud- 
ied here (§ |l|). Ground-based data already revealed 
that CH3OH is not mixed with the bulk of the wa- 



ter ice ( Skinner et al. 1992 ). Recent ISO results indi- 
cate extensive ice segregation toward massive young 
stars; the absorption profiles of CO2 ice indicate that 
the methanol resides in a H2 0:CH3QH:CQ2=1:1:1 layer 
heated to 



80 K as sh own b y [Ehrcnfrcund et al. (1998) 
iGcrakincs et al. (1999)1 and [Dartois et al. (1999a)j The 



¥ 33A 



1 " _ I I □ 

o o o o o & 



GL 2591 



o o olo'fllolo & 



0.1 0.2 0.3 0.4 0.5 



0.1 0.2 0.3 0.4 0.5 



Spatial Frequency (arcsec ) 

Fig. 7. Visibility amplitudes of CH 3 OH 3i -> 4 A+ 
emission as observed with OVRO toward W 33 A (left) 
and GL 2591 (right), and model points for the constant- 
abundance (open circles) and the "jump" models (open 
squares). 
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column density of solid CH3OH is typically ~ 5% of that 
of water ice, but in two cases, W 33A and GL 7009S, it 
is as high as 30%. It has been suggested that CH 3 OH 
ice is preferentially formed around more massive young 
stars, but among the sources studied here, these two do 
not stand out in luminosity. Instead, they are the two 
most embedded sources in the sample (see van der Tak 
et al. 2000), which may enhance the formation of solid 
methanol. The gas-phase abundance of CH 3 OH (Table ||) 
is much lower in GL 7009S than that in W 33A, suggesting 
that GL 7009S is in an earlier evolutionary state where ice 
evaporation affects only a small part of the envelope. 

A link between solid and gaseous methanol is plau- 
sible because the sources with high gas-phase CH 3 OH 
abundances are also the ones with high fractions of an- 
nealed solid 13 C02 (Boogert et al. 2000), and high abun- 
dances of gas-phase H2O and CO2 (Boonman et al. 2000). 
These molecules evaporate at ~ 90 K, like methanol. 
The CH3OH abundances also follow the ratios of enve- 
lope mass to stellar mass and the 45/100 fim colours from 
van der Tak et al. (2000), which confirms the picture that 
warmer sources have higher molecular abundances in the 
gas phase. These results indicate that the excitation and 
abundance of gaseous CH3OH can be used as evolutionary 
indicators during the em bedd ed stage of massive star for- 
mation. As discussed in § 3T, using submillimeter data to 
trace evolution has the advantage of being independent of 
source orientation or total mass, because the dust emission 
is optically thin at these wavelengths. 

Grain mantle evaporation appears to be much less im- 
portant for H2CO than for CH3OH. We have compared 
the H2CO data to the constant-abundance models from 
van der Tak et al. (2000) in similar plots as Figure [|, and 
found good agreement. There is no evidence for jumps in 
the H2CO abundance by factors >3 within the tempera- 
ture range of 20 — 250 K that the observations are sensitive 
to. This result suggests that the formaldehyde observed in 
these sources is predominantly formed in the gas phase by 



In eleven regions of (mostly massive) star formation, 
Mangum & Wootten (1993) derived ratios of ortho- to 



oxidation of CH3. Gas-phase models by [Lee ct al. (1996) 
indicate an abundance of ~ 10~ 9 , similar to the observed 
value, although a contribution from ice evaporation at the 
10~ 9 level cannot be excluded. 

The high observed HDCO/H2CO ratios and the de- 
tections of D 2 CO toward the Compact Ridge in Orion 
(Wright et al. 1996; Turner 1990) and toward embedded 
low-mass objects such as IRAS 16293 (van Dishoeck et 
al. 1995; Ceccarelli et al. 1998) indicate H2CO forma- 
tion on grains. In our sources, HDCO is not detected to 
H 2 CO/HDCO > 10; in a survey of hot core- type sources, 
Hatchcll ct al. ( 1 998 )| obtained DCN/HCN - 10" 3 , much 
lower than in embedded low-mass stars. Although the sur- 
face chemistry should qualitatively be the same, the cold 
(~ 10 K) phase may last too short in our sources to build 
up large amounts of deuterated molecules such as those 
seen in the Compact Ridge and in low-mass objects. 



para-H 2 CO of 1.5 — 3 and took this result as evidence 
for grain surface formation of H2CO. However, this con- 
clusion appears tentative, since for most of their sources, 
only one line of ortho-B^CO was observed. In addition, 
at the H2CO column densities of ~ 10 14 cm -2 , represen- 
tative of the regions studied by Mangum & Wootten and 
by us, the K = and K = 1 lines of H2CO have opti- 
cal depths of ~ 1, requiring careful modeling. The data 
and models presented in this paper are consistent with an 
ortho-/para-H 2 CO ratio of 3. 



High abundances of methanol as observed in the ices, 
~ 10~ 6 relative to hydrogen, cannot be produced in the 
gas phase, except maybe in shocks. This mechanism has 
been pro posed for H2O ice by Bergin et al. (1999) and for 



Charnley fc Kaufman (2000) ; its application to 



C0 2 by 

CH3OB depends on the existence of a high-temperature 
route to form methanol, which is not yet known. 
Hartquist et al. (1995) proposed the reaction CB4+OH, 



but ISO-SWS observations by [Boogert et al. (1998)| indi- 
cate low abundances of gaseous methane in our sources, 
and at T>200 K, all OB should be consumed by B and 
B 2 to form B 2 0. The widths of the CB 3 OB lines of only 
a few km s _1 also argue against formation in shocks. 

More likely, the solid methanol is formed by reactions 
on or inside the ice layers. Addition of B atoms to CO 
molecules will lead to B2CO, and further to CB3OB. In 
the literature, there are three proposed sources of atoms: 
direct accretion from the gas phase, ultraviolet irradiation 
and bom bardment by energetic particles. We show in Sec- 
tion 5A that the latter two mechanisms are unlikely to be 
important for the sources studied in this paper, and focus 
for now on the first. 



5.2. Surface chemistry model 

If the CB3OB ice that we see evaporating in these sources 
originates from B atom addition (=reduction) of CO ice 
on the surfaces of dust grains, it cannot have been pro- 
duced under the current physical conditions. The evapo- 
ration temperature of CB3OB and B2O ice, 90 K, is much 
higher than that of CO and O, w 20 K, and that of H, 
~ 10 — 15 K, depending on its surface mobility and reactiv- 
ity. Any formation of CB3OB ice through surface reduc- 
tion of CO must therefore have occurred before the central 
star heated up its envelope above ~ 15 K. The most likely 
phase of the cloud to form methanol ice through surface 
chemistry is therefore the pre-protostellar phase, when the 
grain temperature may have been as low as ~ 10 K and 
the cloud was contracting to form a dense core. Bowever, 
about 1/3 of solid CO is observed inside the water ice 
layer, and will not evaporate until mu ch higher tempera- 
tures are reached ( Tielens et al. 1991 ) . If B and O atoms 
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by the curve in Fig. ||, taken from |Charnlcy ct al. (1997) 
with the reaction probability ratio </>h of CO + H to H2CO 
+ H is 1CP 3 . Along the curve, the H/CO abundance ratio 
decreases down to 10~ 2 at its tip. At the density in the 
central regions of our sources, ~ 10 6 cm -3 , atomic H is 
mainly produced by cosmic-ray ionization of H2 , giving a 
constant concentration of ~ 1 cm~ 3 , or an abundance of 
10 -6 . The abundance of CO in these sources was measured 
by van der Tak et al. (2000) to be « 1 x 10~ 4 , as expected 
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if all gas-phase carbon is in CO and no significant deple- 
tion occurs. Hence, H/CO « 10 -2 , and our sources should 
lie right at the tip of the curve. The value </>h = 10 -3 
implies that the CO^ H2CO reaction limits the rate of 
CH3OH formation, consistent with the non-detection of 
evaporated H 2 CO. 

The figure also plots the observed abundance ratios, 
where the H2CO and CO abundances have been taken 
from van der Tak et al. (2000), while for CH3OH, either 



the constant-abundance model or the "top level" of the 
jump model has been plotted. The data points are seen to 
lie along a "sequence" which lies along the model curve, 
perhaps displaced toward lower CH3OH abundances by an 
order of magnitude. Although this behaviour is consistent 
with solid-state formation of H 2 CO and CH 3 OH at vari- 
ous H/CO ratios, the different rotation temperatures and 
abundance profiles of H2CO and CH3OH argue against a 
common chemical origin for these two species, and favour 
complete hydrogenation of CO into CH3OH. The figure 



CO/CH3OH 

Fig. 8. Filled circles: observed abundance ratios of 
H2CO/CH3OH (top) and CO2/CH3OH (bottom) versus 
that of CO/CH3OH. Open triangle: data on translucent 
clouds from [Turner (1998)]; open square: data on the com- 
pact ridge in Orion from Sutton ct al. (1995). Stars: solid- 



state data from Keane et al. (2000) 



Heavy line: surface 
(1997)| . Open cir- 



chemistry model after Charnley et al 
cles: results of surface chemistry model 
and density of the gas phase. 



can be similarly trapped, solid-state reduction and oxida- 
tion may occur at temperatures well above 20 K. 

Tielens & Hagen (1982) and Tielens & Allaman- 
dola (1987) proposed that direct accretion of atoms and 
molecules from the gas phase, followed by low-temperature 
surface reactions, determine the composition of grain man- 
tles. Activation barriers are offset by the long effective du- 
ration of the collision on the surface. The resulting mantle 
composition is determined by the relative accretion rates 
of H and CO onto the grains, and the relative height of 
the reaction barriers. 

This mechanism predicts that the abundance ratios of 
H2CO/CH3OH and CO/CH3OH are correlated, as shown 



also plots solid state observations by Keane et al. (2000) 
which are lower limits in the case of CO because of evap- 
oration. 

Since the abundances of CO and H 2 CO are almost the 
same in all our sources, the spread in the abundance ra- 
tios is probably due to different degrees of CH3OH evapo- 
ration and subsequent gas-phase destruction. Indeed, the 
abundances of gaseous CH3OH are factors of > 10 lower 
than the values measured in the solid state. This differ- 
labeled by time ence is expected since methanol is destroyed in ~ 10 4 yr 
(Charnley et al. 1992) in reactions with ions such H3", 
leading to a rich chemistry with species like CH3OCH3 
and CH3OCHO. These species are indeed detected toward 
the "hot cores" W 3 (H 2 0) and NGC 6334 IRS1, weaker 
toward the "jump sources" GL 2591 and NGC 7538 IRS1, 
and absent in the low-methanol source W 3 IRS5 (Helmich 
& van Dishoeck 1997; van Dishoeck & van der Tak 2000). 

We conclude that the observed spread in the gaseous 
methanol abundance is due to incomplete evaporation 
for the cold sources, which are the least evolved ones in 
our sample, and to gas-phase reactions for the hot cores, 
which are the most evolved. The most likely source of the 
methanol is grain surface chemistry in the pre-protostellar 
phase. The conditions in this phase cannot be derived by 
comparison to CO and H 2 CO, for which our data suggest 
that gas-phase processes control the chemistry. A more 
promising molecule is CO2. 
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5. 3. The H/ O competition: Density and duration of the 
pre-stellar phase 

The ISO detection of ubiquitous solid CO2 by 



Table 4. Results of gas-grain chemical model at T — 10 K. 



Gerakines et al. (1999) makes clear that oxidation of CO 
is a potentially important competitor to reduction. These 
data are shown in the bottom pa nel of Fig. p|, as well as 
gas-phase CO2 abundances from Boonman et al. (2000) 



Gas-grain chemistry is clearly important for CO2, as it 
is for CH3OH, warr anting a comparative study of the 
two. The model by Charnley et al. (1997) did not in 
elude reactions with O, while Tielens fc Hagen (1982 



used O2 as the dominant form of oxyge n, which is in - 
consistent with recent observational limits ( Melnick 2000|) . 
We have constructed a model of gas-grain chemistry 
based on the modified rate equation approach described 
in Casclli et al. (1998) , and extensively tested against a 
Monte Carlo program. The chemical system consists of 
three gaseous species: H, O and CO, and eight surface 
species: OH, H 2 0, 2 , C0 2 , HCO, H 2 CO, CH 2 OH and 
CH3OH. The calculations assume fast quantum tunnel- 



ing of hydrogen, as proposed by Tielens & Hagen (1982) 
If hydrogen docs not tunnel, as measurements by 
Katz et al. (1999) suggest, our simple system is unable 
to form significant mole fractions of CH3OH and CO2. 
In the case of reduced mobility for all species, the Eley- 
Rideal mechanism should be considered to model surface 



chemistry (Hcrbst 2000), which is not included here. 

The mantle composition (Table |]) depends on den- 
sity through the composition of the gas phase, in partic- 
ular the H/O ratio. We take t he concentrations of H, O 
and CO from Lee et al. (1996) , for the case of low metals, 
T = 10 K, in steady state. The main result is that after re- 
duction to HCO, CO is mostly reduced into CH3OH at low 
densities and mostly oxidized at higher densities. In our 
model, the reaction CO+O— >C02 has a barrier of 1000 K 



(see d'Hendecourt et al. 1985), else ices of 60% CO2 would 
form which are not observed. Water is copiously made 
through the sequence O+H^OH and OH+H^H 2 0. The 
detailed composition of the grain mantles is quite uncer- 
tain because of unknown reaction rates. For instance, the 
fraction of H2CO depends on relative barrier height of the 
CO+H and H2CO+H reactions, but H2CO is never a ma- 
jor component of the ice layer, consistent with gas-phase 
and solid state observations. The results in Table || refer 
to a chemical time scale of 1000 years. Unlike inert species 
such as O2, H2O, CH3OH and CO2, the fractions of trace 
species (fractions < 0.01) also depend somewhat on time, 
and are not used in the analysis. Species that react quickly 
with H maintain a constant population because they are 
hydrogenated on the H accretion time scale (~ 1 day if 
ns = 10 3 cm -3 ), and their molar fractions decline as the 
ice layer grows with time. 

As an example, we will now estimate the yield of 



CH3OH by this scheme for the case nu 



10 4 



Species Total density tih 



1(T 10 4 10 s 



Gas phase concentrations (cm 3 ) 



H 


1.15 


1.15 


1.1 


O 


0.09 


0.75 


7.0 


CO 


0.075 


0.75 


7.5 


Molar fractions 


in the solid state 


2 


9.3(-3) 


0.11 


0.27 


H 2 


0.60 


0.37 


0.12 


H 2 


0.0 


0.0 


0.0 


CO 


1.4(-5) 


3.9(-5) 


0.48 


H 2 CO 


1.4(-5) 


3.1(-5) 


3.1(-2) 


CH3OH 


0.39 


0.41 


7.8(-3) 


C0 2 


6.1(-3) 


0.12 


8.7(-2) 



is fd — ndcVcoS, with nd and a the number density and 
cross section of the dust grains, Vco the thermal velocity 
of CO and S the sticking coefficient, assumed to be 10%. 
Taking silicate grains of radius 0.1 fim, fd = 3 x 10~ 15 s _1 , 
or 10% in 10 6 years. At this density, 77% of this de- 
pleted CO goes into CH3OH, so that upon evaporation, 
the abundance ratio would be CH3OH/CO w 0.086 or 
CH 3 OH/H 2 - 10~ 5 . 

Assuming that the composition of ice layers is deter- 
mined by surface chemistry, our model can be used to 
investigate the initial conditions of massive star forma- 
tion by considering abundance ratios of CO, C0 2 and 
CH3OH. The CO2/CH3OH ratio is sensitive to density, 
independent of time, and equals the ratio of the molar 
fractions of solid C0 2 and CH 3 OH in Table @. Time is 
measured by the ratio of raw material (CO) to product 
(CH3OH) through the fraction of CO in the solid state de- 
rived above, 0.01(i/10 5 yr)(nH/10 4 cm -3 ). The abundance 
ratio in the gas phase is this fraction multiplied by the 
efficiency of CO^CHsOH conversion, equal to the mo- 
lar fraction of CH3OH divided by the sum of the frac- 
tions of CO and its possible products, CO2, H 2 CO and 
CH3OH. Figure H plots the synthetic abundance ratios as 
open circles, labeled by time and density. If this model 
is valid, the observed CO2/CH3OH ratio constrains the 
density in the pre-protostellar phase to be w 10 5 cm~ 3 , 
or higher, since some CO2 may be destroyed shortly af- 
ter evaporation. At this density, hydrogenation of CO is 
incomplete due to the low H flux, and a significant abun- 
dance of solid H2CO is expected, consistent with the re- 
sults of Keane et al. (2000)| . Using this limit on the den- 



The depletion rate of CO molecules from the gas phase 



sity, the observed CO/CH3OH ratio constrains the time 
spent in the pre-protostellar phase to be < 10 5 years within 
the grain surface chemistry scenario. The same conclu- 
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sion is reached when the CO/CO2 ratio is used instead 
of CO/CH3OH; the uncertainty is a factor of ten due 
to the unknown sticking coefficient. This time scale is 
significantly smaller than the corresponding number for 
low-mass stars, where this phase lasts ~ 10 6 years as de- 
rived from the ratio of dense cores with and without stars 



found by Bcichman et al. (1986). The same time scale of 
~ 10 6 yr is obtained from our models for an abundance 
ratio of CO/CH3OH ~ 10 3 , as observed in the compact 
ridge in Orion. 

These models also help to understand observations of 
solid CO2 in other regions where no star formation is oc- 
curring. Toward the field star Elias 16 behind the Taurus 



dark cloud, CO2 has been detected (Whittct et al. 1998 ) 
but CH3OH has not flChiar et al. 1998|) . The models sug- 
gest thus that the density in this region is > 3 x 10 4 cnr 3 . 
The model does not explain the lack o f solid CH3OH to - 
ward SgrA* ( < 3% relative to H 2 ice; |Chiar et al. 2000| ), 
where CO2 ice has been detected. For this low-density 
line of sight, the opposite ratio would be expected. An 
enhanced temperature or ultraviolet irradiation may be 
important here. For the disk aroun d N GC 7538 IRS1, 
where iV(H 2 CO)/iV(CH30H) >1 (§ the models in- 
dicate that the density is > 10 5 cm -3 . 

The model is also consistent with current observational 
limits on solid 2 . At low densities, our models drive all 
O into H 2 on the dust grains, but for > 10 5 cm~ 3 , part 
of it goes into 2 . However, the ratio 2 /CO remains 
< 1 on the surface, in agreement with the limit derived 
by Vandenbussche et al. (1999) from ISO-SWS observa- 
tions of NGC 7538 IRS9. Evaporation of solid 2 plays 
a role, but observations of gaseous 2 with the Submil- 
limeter Wave Astronomy Satellite give an upper limit to- 
ward massive star-forming regions of ~ 10~ 7 relative to 
H 2 ( Mclnick 2000| ), consistent with the solid-state results. 



5.4- Alternative models 

Could ultraviolet irradiation or energetic particle bom- 
bardment also produce the observed trends in the abun- 
dances? The production rate of species i by irradiation 
can be written as drii/dt = a^Airr^ with $ the ul- 
traviolet flux in photons cm~ 2 s^ 1 and r g the radius 
of the dust grains, taken to be 10~ 5 cm. Inside dense 
clouds, most ultraviolet radiation is produced by cosmic- 
ray interaction with H 2 , which gives a field approxi- 
mately equal to the interstellar radiation field at Ay = 5 
( [Prasad fc Tarafdar 1983| ), or 5000 photons cm~ 2 s" 1 . 
The reaction yields on follow from laboratory experiments 
such as those reported at the Leiden Observatory Labo- 
ratory database^, in this case on a mixture of initial com- 
position H 2 0:CO = 100:33. The amounts of CH 3 OH and 
C0 2 produced do not depend on this ratio as long as it 
is > 1. Using band strengths by Gerakines et al. (1996) 



10 8 s 1 per grain, corn- 



see www.strw.leidenuniv.nl/~lab 



and [Kerkhof ct al. (1999)| , we obtain a = 1.9 X 10 3 for 
CH3OH and a = 3.2 x 10~ 2 for C0 2 . The ratio of these 
numbers is in good agreement with the observed abun- 
dance ratio of ~ 10. However, the absolute values of the 
on imply production rates of 
pared to the accretion rate of H, O and CO of ' 
Stellar ultraviolet radiation can only affect the inner parts 
of the envelopes since their extinctions are Ay ~ 100 mag- 
nitudes. The absence of a strong ultraviolet field through- 
out the envelopes of our sources is also suggested by the 
observational limits from ISO-SWS on mid-infrared fine 
structure lines and of emission by polycyclic aromatic hy- 
drocarbons ( van Dishoeck fe van dcr Tak| 2000). 

The processing of interstellar ice by cosmic rays 
was studied by Hudson fe Moore (1999)| , who bombarded 
an H 2 0:CO=5:l ice mixture with protons of energy 
0.8 MeV. These experiments produced abundance ratios of 
H 2 CO/CH 3 OH=0.6 and CO 2 /CH 3 OH=2.0, in fair agree- 
ment with the observed values. Other species formed as 
well, notably formic acid (HCOOH) and methane (CH4), 
which indeed are observed in interstellar grain mantles. 
However, the experiments by Hudson & Moore produced 
almost twice as much HCOOH and CH4 as C0 2 , while 
t he observed abund ances are only rs 10% of that of C0 2 
flgchutte et al. 1999| ; Boogert et al. 1998). In addition, the 
particle dose in the experiments was 2 x 10 15 cm -2 , while 
the interstellar cosmic-ray flux is only 3 cm~ 2 s^ 1 . Hence, 
the experiments simulate a bombardment for ~ 3 x 10 7 yr, 
a factor of ~ 1000 longer than the ages of the sources stud- 
ied here. Stellar X-ray emission (Glassgold et al. 2000) 
only acts on small scales, especially because of heavy at- 
tenuation in the envelopes of these obscured objects. 



6. Conclusions 

The chemistry of CH3OH and H 2 CO in thirteen regions 
of massive star formation is studied through single-dish 
(JCMT) and interferometer (OVRO) line observations at 
submillimcter wavelengths. Our main conclusions are: 

1. The submillimcter emission lines of CH3OH toward 
most sources have widths of 3 — 5 km s _1 , consistent with 
those found earlier for C 17 and C 34 S. However, in the 
low-luminosity sources GL 7009S and IRAS 20126, the 
line shapes reveal that CH3OH is present in the outflow. 
These results indicate that the desorption of ices in the 
envelopes of low-mass protostars is primarily by shocks, 
while thermal processes dominate in the case of massive 
stars. 

2. Rotational temperatures of CH3OH range from 10 
to 200 K and correlate very well with the excitation tem- 
perature of C 2 H 2 measured in infrared absorption. This 
correlation suggests that both species trace the same gas, 
which may be outgassing of icy grain mantles. For H 2 CO, 
the range in T rot is only 60 — 90 K without relation to 
?ex(C 2 H 2 ), suggesting a different chemical origin. We pro- 
pose that T rot (CH3OH) can be used as evolutionary indi- 
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cator during the embedded phase of massive star forma- 
tion, independent of source optical depth or orientation. 

3. Detailed non-LTE radiative transfer models of 
the CH3OH lines suggest a distinction of three types 
of sources: those with CH 3 OH/H 2 ^ 10~ 9 , those with 
CH 3 OH/H 2 ^ 1CT 7 and those which require a jump in 
its abundance from ~ 10~ 9 to ~ 10~ 7 . The models are 
consistent with the w 3" size of the CH 3 OH 107 GHz 
emission measured interferometrically. The location of the 
jump at T « 100 K strongly suggests that the methanol 
enhancement is due to evaporation of icy grain mantles. 
The sequence of low-methanol — > jump — > high-methanol 
sources corresponds to a progression in the ratio of enve- 
lope mass to stellar mass and the mean temperature of 
the envelope. In contrast, the observed H 2 CO seems pri- 
marily produced in the gas phase, since the H2CO data 
can be well fit with a constant abundance of a few x 10~ 9 
throughout the envelope. The grain surface hydrogenation 
of CO thus appears to be completed into CH3OH, with lit- 
tle H 2 CO left over. 

4. Model calculations of gas-grain chemistry show 
that CO is primarily reduced (into CH3OH) at densi- 
ties 71h~ 10 cm -3 , and primarily oxydized (into C0 2 ) at 
higher densities. To keep sufficient CO on the grains, this 
mechanism requires temperatures of < 15 K, i.e., condi- 
tions before star formation. Assuming that surface reac- 
tions proceed at the accretion rate of CO, the observed 
C0 2 and CH3OH abundances constrain the density in the 
pre-protostellar phase to be tih ~ a few 10 4 cm~ 3 , and the 
time spent in this phase to be < 10 5 yr. Our surface chem- 
istry model predicts that lines of sight through clouds with 
a high H/O ratio will show abundant solid methanol and 
less C0 2 . Ultraviolet photolysis and radiolysis by ener- 
getic (MeV) protons appear less efficient as ice processing 
mechanisms for these sources; radiolysis also overproduces 
HCOOH and CH 4 . 
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